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Abstract 

Quadratic systems generated using Yang-Baxter equations are integrable in a sense, 
but we display a deterioration in the possession of the Painleve property as the number 
of equations in each 'integrable system' increases. Certain intermediate systems are 
constructed and also tested for the Painleve property. The Lie symmetries are also 
computed for completeness. 

1 Introduction 

Golubchik and Sokolov [17] recently discussed systems of ordinary differential equations 



which are integrable by the standard factorisation method of Adler-Kostant-Symes |2l[] 
or the generalised factorisation method fT6|| . They established relationships between such 
reductions, operator Yang-Baxter equations and some kinds of nonassociative algebras. In 



their paper [17] a number of specific examples was given without much detailed comment 
on the solutions except for the possibility of the system passing or not passing the Painleve 
test. The examples are as follows: 

P t = P 2 -RP- QS 

Qt = (f3- 2)RQ + (3PQ 

Rt = R 2 - RP -QS 

S t = {3- [5)RS + (1 - 0)PS 



Pt = RQ 

Qt = PQ (1.2) 
R t = aRP 
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P t = 2PR + XQR 

Qt = 2QR - XPR (1.3) 
R t = P 2 + Q 2 + R 2 

P t = (y- n)QR 

Qt = 2fiPQ + iiQ 2 + vQR (1.4) 
R t = -2vPR - vR 2 - fiQR 

— = X 2 + 3XY = X(X + 3Y) 

— = 3XY + Y 2 = Y(3X + Y) 

in which Greek letters are constants. We have given the sets of equations in the order and 



and 



the form as presented by |17]. 



The integr ability of these systems, apart from the two-dimensional system in (jl,5|), is 
not immediately obvious, although in the paper of Golubchik and Sokolov |o] we are 
informed that they are so constructed to be integrable in the sense of the generalised 
factorisation method. Unfortunately they do not connect this definition of integrability 
with the usual ones which are as follows: 

(a) The existence of an explicit solution relating the variables, not necessarily the depen- 

dent variables as explicit functions of the independent variable, or, in the instance 
of systems, a number of explicit independent functions equal to the number of the 
dependent variables; 

(b) The existence of a sufficient number of independent first integrals and invariants 
which could be used to give a local version via the implicit function theorem of (a) 
above; 

(c) The existence of a sufficient number of symmetries either to reduce the differential 

equations of the system to algebraic equations or to obtain the independent first 
integrals and invariants of (b) above; 

(d) The passing of the Painleve test, either in the normal form or the weak form, which, in 

the form as given by Conte || , will guarantee the existence of a solution as required 
by (a) above that is analytic in the complex t-plane or in a region of it apart from 
isolated moveable singularities which are poles or algebraic branch points. 

The only hint which the two authors give is their comment associated with fll.lD to the 
effect that for general values of the parameter (3, one would not expect the system to 
pass the Painleve-Kovalevskaya test. One infers from this that they are not looking for 
a solution which is analytic apart from isolated singularities which are either poles or 
algebraic and so one can presume that they intend (a) above. 

There are those who would insist that a solution be a (single-valued) function (cf Conte 
Hi) and, in the strict mathematical sense, they are quite correct. However, in practice the 
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domain of the variable, t, in which the solution is going to be used is much smaller than 
the complex t-plane. It is with this emphasis on the practicality of applications that (a) 
above is proposed. 

It is the intention of this paper to take up the systems enumerated above where the 
authors of [17] left them. We shall see that the solutions of these systems of equations have 
certain interesting aspects which illustrates various properties of differential equations. 

The most obvious one is that all of the systems quoted are quadratic systems which 
are of tremendous interest in many areas of mathematical modelling ranging from studies 
of population systems to the reactions in chemical systems. They, along with the closely 
related Lotka-Volterra systems, are usually studied from the point of view of dynamical 
systems since in general they are not integrable, but are well adapted for intensive and 
informative qualitative investigation. Nevertheless there has been a number of studies 
of these systems and the circumstances under which they possess first integrals and/or 
invariants ||, 0. The connection between the possession of the Painleve property and 
the existence of first integrals/invariants has also been investigated (jOJ. 

A less obvious property of the systems presented by Golubchik and Sokolov is that 
all of them possess the two point symmetries of invariance under time translation and 
the selfsimilar transformation. These two symmetries immediately make them interesting 
from the point of view of the Painleve property since there is a very close relationship 
between the nature of the singularity and the specific form of the selfsimilar symmetry 
[13]. The selfsimilar symmetry permits an analysis of the next to leading order behaviour 
and the Painleve property follows from this when the next to leading order behaviour 
becomes a Laurent expansion about the singularity [14]. 

There is the question whether the possession of the two symmetries mentioned above is 
sufficient to guarantee integrability. As we have been informed in [17], all of these systems 
are integrable, we have an excellent opportunity to examine the symmetry structure of 
these systems and equations related to them in the knowledge that the results will char- 
acterise integrable systems. Of course it is an open question whether the examination 
will lead to a clear result. We defer our comments on the question until the conclusion. 
However, our path to the conclusion will reveal a number of interesting features of these 
systems and their connections to other well known facets of integrable systems. 

On the presumption that complexity of appearance is a representative of complexity of 
fact we shall more or less work backwards through the list of systems above and attempt 
to examine them in increasing order of difficulty. In the case of some of the systems we 
make use of rescaling to make them look simpler. This may simply be a psychological 
trick, but it is one worth practising on all occasions. One of our principal activities will 
be the examination of these systems and any derivate equations for their Lie symmetries. 
As an aid in this examination we make use of Program Lie [18, 33 1 to attend to the dirty 
business of the calculation of the symmetries. This program, which is one of the more 
successful analysers for symmetry, is freely available on Simtel sites. Naturally we shall 
also be concerned with the possession of the Painleve-Kovalevskaya property either by a 
direct computation or by a consideration of the structure of the solution of the system of 
equations. 

Before we commence our analysis we mention some terminology introduced in reference 
[24]. In systems of differential equations invariant under time translation and rescaling 
all terms contribute to the dominant behaviour. However, it is possible that a subset of 
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terms has an additional or different rescaling property. This symmetry is termed 'sub- 
self similarity' and the terms possessing this symmetry are called 'sub dominant'. These 
terms also have to be analysed for possession of the Painleve property, naturally with 
compability at the resonances no longer guaranteed. A simple example of the possession 
of subselfsimilarity is the generalised Chazy equation 

x + axx + bx 2 = 0. (1-6) 

The equation as a whole has the two symmetries G\ = d/dt and G2 = —td/dt + xd/dx. 
The second and third terms have the two homogeneity symmetries td/dt and xd/dx and 
so would also have to be analysed separately. 

Finally we note that this is the first in a series of papers devoted to the intertwined 



subjects of symmetry, singularity and integrability. The second paper in the series [25] 
treats the singularity analysis and integrability of the class of second order equations 
invariant under time-translation and rescaling and the third [26] some aspects of complete 



symmetry groups |22J. 



2 The two-dimensional system 

2.1 Painleve analysis 

We rewrite (|l.5|) , viz 

^ = X 2 + 3XY 

dt (2.1) 
dY 

= 3XY + Y 2 

at 



MS 



x = x 2 + 3xy (2.2a) 
y = 3xy + y 2 . (2.2b) 



When we perform the Painleve analysis on the system fl2.2| ), we find that, using the 
usual notations and methods, the singularity corresponding to the selfsimilar symmetry 
of the system gives the following behaviour 

x = -\ T ~ l + 3/ir"2, p = -l, r = -l,| 
V = ~\r~ X ~ 2a*t-2, q = -l, r = -l,i, 

where r = t — to an d to is the location of the movable singularity and fi is the second 
arbitrary constant of integration, which indicates the weak Painleve property implied in 
the solution (see ( 2.31| )). 



There is also the possibility of a different type of singular behaviour for ( p.2|) which 
would be given by the two equivalent possibilities p = — l,q > — 1 and p > —l,q = — 1. 
This type of behaviour is found for example in the analysis of the Mixmaster universe 
1C[. A simple-minded approach would be to write 

00 00 

X>r-\ y = 5>* (2.4) 

i=0 i=0 
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(2.5) 



for the first of these possibilities. This yields the particular solution 
-1 

x = — , y = 0. 

T 

However, bearing in mind the expansion in powers of r2 found in ( |2.3[) we do well to 
propose 

oo , oo 



X 



E' 

i=0 



i-l 



y 



E 

i=0 



i . i 



T2' 



(2.6) 



This gives the one-parameter solution 



13 
54 
1 



+ ... 



(2.7) 



Since ( |2.7| ) is a one-parameter system, this pattern of singular behaviour (and clearly the 
one corresponding to the interchange of the roles of x and y) does not pass the weak 
Painleve test and so the system (2.2) does not possess the weak Painleve Property. (In 
this we follow the definition of Tabor [p4| [p 330].) 



2.2 First reduction 



We now obtain the solution of (2.2) as follows. Firstly we are immediately struck by the 
resemblance to the coefficients in the expansion of a cubic binomial. This suggests the 
change of variables 



2 , 2 

u = x + y 



v = x + y 



so that the system (2^2) takes the form 



u = 2v 6 



-2u + 3v z 



(2.8) 

(2.9a) 
(2.9b) 



which by differentiation of the second of these and substitution for u gives the second 
order equation 



6vv + 4u d = 0. 



(2.10) 



Equation fl2.10| ) is of a type [20, 6.43, p 551] which has been studied extensively in the 
literature because of the frequency of its occurrence in different areas Jj], |l2], [l5|, |3l], 32] 
and in particular is known, for certain values of the coefficients, to have eight symmetries 
instead of the obvious two of invariance under time translation and selfsimilarity, thereby 



making it linearisable |3C]. 

For the system (2.9) we obtain two possible sets of leading order behaviour. They are 
p = —2 and q = —1 with the resonances r = —1,1 and p = q = — 1 with the resonances 
r = —1,-2. We give the critical parts of the expansions. They are 



It- 2 +3HT- 1 + ... 
i T -l_ 2/i + ... , 



(2.11) 
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which gives the expansion about the singularity at t = to, and 

(2.12) 



u = t~ 1 + 2fit~ 2 + 3vt~ 3 



v = t" 1 + pdT 2 + vt 3 , 
which represents the asymptotic expansion of the solution. The second constant of inte- 



gration is v and we have followed the practice introduced by Feix et al [13] of making the 
expansion in the variable t because of the asymptotic nature of the expansion. 

The expansion in ( 2.11| ) has been termed [|lj| a right Painleve series because it is a 



Laurent expansion about the singularity at to i n ascending powers of r whereas that in 
( 2.12] ) is called a left Painleve series because it is a Laurent expansion in descending powers 



of the independent variable, t. The existence of what is now called the left Painleve series, 
that is an asymptotic series as a representation of the solution, has only been recognised 
in recent years j27j. 

It is interesting that the subselfsimilarity symmetry property of the system (|2.2| ) and 
its corresponding Painleve analysis disappears under the transformation which leads to 
the system ( |2.9D . The reason for this is obvious for the first of (|2.9| ) contains only two 
terms and a minimum of three terms is required for the presence of a subselfsimilarity 
symmetry [|L| (these can be additive or multiplicative). It is an interesting speculation 
as to whether there is always a transformation which will remove the possibility of the 
subself similar behaviour. 



2.3 Second reduction 



An alternative route to the analysis of (2^2) initiated by the transformation ( |2i 
by the transformation 



x + y, 



v = x 



The transformed system (|2.2|) is now 



u = 2u z 
V = uv. 



is given 



(2.13) 



(2.14a) 
(2.14b) 



We have two possible routes for the solution of ( 2.14 ), viz the differentiation of ( 2.14 a) 
and elimination of v or the differentiation of ( 2.14| b) and elimination of u. The first route 
produces the nonlinear second order ordinary differential equation 



2uii - Avr 



(2.15) 



which belongs to the same class of equations as ( [2.10 ). If we perform the leading order 
analysis on ( 2.15D , we find that the only possible singularity is a simple pole since there 
is no possibility of subselfsimilarity. The coefficient of the leading power can be either 
— 2 or 1. In the former case the resonances are at r = —1,3 and in the latter case at 
r = —1,6. There is no possibility of inconsistency at the second resonance because of the 
rescaling symmetry. Both expansions are right Painleve series and the equation possesses 
the Painleve property. If we look to the solution of ( 2.15j ) by means of reduction of order, 
we come to an Abel's equation of the second kind and firstly transformation is necessary 
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to obtain the solution. An alternative route is to use the Riccati transformation of ( [2.27 ) 
(with the same value of a) to obtain the third order equation 



w w 



2ww = 



(2.16) 



which has a solution in terms of elliptic functions. 

The second route produces the nonlinear second order ordinary differential equation 



vv 



3v 2 + v 4 







(2.17) 



which belongs to a class of equations to be found in Kamke [20, 6,128, p 574]. Under the 



transformation v = w 2 ( 2.17j ) takes the particularly simple form 



w 



(2.18) 



which we shall again encounter in (3.5). The solution of ( 2.18; ) is trivial and so we have 

i 



v(t) = (A + 2Bt + t z y2. 



(2.19) 



The other solution u(i) follows easily from ( [2.14| b). One easily checks that these two 
solutions lead to the same forms for the general solutions of the original system, ( |2.2[ ). 

Equation ( p. 15 ) has just the two Lie point symmetries of invariance under time trans- 
lation and selfsimilarity. In marked contrast ( [2.17 ) possesses eight Lie point symmetries 
which, in a sense, makes it more like ( |2.1C ) than (2.15) is. We observe that the first two 
terms of ( 2.17 ) are invariant under the homogeneity symmetry ud/du and so expect a 
more complicated pattern of leading order behaviour. This is borne out by the analysis 
which gives in the case of all terms dominant an exponent of —1 with coefficients ±1 and 
in the subdominant case exponents of and — ^. In the former case the resonances are 
given by r = —1, —2 and so for this type of leading order behaviour the equation passes 
the Painleve test and has a left Painleve series, that is, an asymptotic solution. This is 
easily seen from ( |2.19| ) if we expand it in the asymptotic form 



/ N 1 ( B 

v (t) = - 1 

w t \ 2t 



W 2 - 4A 
8*2 



(2.20) 



For the subdominant behaviour and p = the first few terms of the Taylor expansion give 



v(t) = a + ait + 



2a 



(3a 2 



at) t 2 + 



1 — 3a, 



ait 3 + 



(2.21) 



We observe that the third term in ( 2.17 ), which breaks the homogeneity symmetry of the 
first two terms, does not affect the existence of the Taylor series solution. The coefficients 
of the first two terms are arbitrary constants and so this is the series representation of a 
general solution. In the case that p = 



we obtain 



v (t) = a r 



1 3 5 



5„4 



7 



a^r22 + ... . 



8"0 



(2.22) 



This series representation of the solution contains the two arbitrary constants ao and to- 
The equation satisfies the requirements of the weak Painleve test. This is one occasion 
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when all possible leading order behaviours possess the Painleve property and so there 
is no question that the equation ( 2. 17f) possesses the Painleve property. For the sake of 



completeness, not to mention a little amusement, we perform the Painleve analysis of the 



system ( 2.14 ). With the usual substitution for the leading order behaviour the exponents 
are 

P ~] 2P 2q (2.23) 
q- 1 p + q V ' 

so that we have the possibility of all terms being dominant with p = q = — 1 or of just one 
variable being dominant with p = — l,q > — 1. In the former case the resonances occur at 
r = —1, —2 and so we have a left Painleve series. Because of the symmetry of the system 
we shall not have inconsistencies at the resonances and we obtain the first few terms of 
the expansion as 

U \ , f 1 "\ ...-2 , ( ~W \ ..,-3 



vj-y p ) l + {-(3 )>« + { i ( 2 - 24 ) 

where \i and v are the two arbitrary constants required for a general solution. For the 
second case the substitution of the series expansion 

oo oo 
i=0 i=0 

results in the particular solution 

u(t) = -— v(t) = 0. (2.26) 
2r 



We thus arrive at the interesting conclusion that, while we started with a system, (2.2), 
which does not possess the (weak) Painleve property, we can map it through the nonlinear 
transformation (|2.8|) to the new system (|2,9| ) that has the strong Painleve property and 



through the linear nonsingular transformation ( 2.13| ) to the new system ( |2.14| ) which 



although it has the (strong) Painleve property, naturally splits into two single equations, 
( 2.15| ) and ( 2.17| ), one possessing the (strong) Painleve property and the other the weak 



viz 



Painleve property. 



2.4 General solution 

It so happens that ( [2.10D has the combination of coefficients which gives eight Lie point 
symmetries. In principle one can use these symmetries to obtain the solution, but they are 
particularly complicated and there is an easier route to the solution. It is known that 
the equations which belong to this class can be expressed as a simple third order equation 
by means of a Riccati transformation. We put 



w 

v = a—, 

w 



(2.27) 
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where a is a constant to be determined to give the third order equation an optimal 
simplicity [|J. Using ( 2.27|) we find that ( 2.10| ) becomes 



w , . ww , 9 , w 3 

- - 3 + 6a — + 2 + 6a + 4a 2 = 

u; w z w 6 (2.20) 

45 w = 



when we put a = — \ for both of the terms containing a then vanish. Equation ( 2.28 ) is 
trivial to solve and we obtain 

B + Ct 

"«> = - A + 2Bt + CV <2 - 29) 



where A, B and C are constants of integration. Equation ( |2.9| a) is now a simple quadrature 
and we obtain 

_ gcg + ct) 2 - (g - ac) 

UW " (A + 2 J Bt + Ct 2 ) 2 ' 1 ^ 

for which an additive constant of integration must be set at zero to maintain consistency 
with ( |2.9| b). (Despite the presence of the three constants A, B and C they provide only 
two independent constants of integration.) Finally we obtain the solutions of the original 
system of equations. We have 



l B + Ct±^{B + CtY-2{B^-AC) 
X[t) ~ 5 A + 2Bt + Ct 2 [ 6 ) 

and y{t) has the same form with the opposite sign attached to the square root. We observe 
that in general the solutions x(t) and y{t) have branch point singularities whereas the only 
singularities of u(t) and v(t) are always poles. 

2.5 Symmetry analysis 

We complete our considerations of the system (|2.2| ) by looking at its first integral and 
invariant and some of its Lie point symmetries. Since the systems (|2.2|) and fl2.9|) (or 



( |2.2| ) and ( 2.14 )) are related by a point transformation, the discussion of one system is 
equivalent to the discussion of the other system when it comes to Lie point symmetries 
and, as these underlie the existence of the first integral and invariant, of them as well. 
We have noted that all of the systems discussed in this paper have the two Lie point 
symmetries which leave them invariant under the transformations of time translation and 
rescaling, that is 

Gl = | and = (2.32, 

where the summation over i includes all of the dependent variables. The first integral is 
associated with Gi and is simply obtained from the integration of the ratio of (|2.2| a) and 
(IJb), viz 

dx x 2 + 3xy 

— = %. 2.33 

ay 3xy + y £ 
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It is 



(x - yf 



xy 



(2.34) 



We use G2 to obtain the invariant from the autonomous first order ordinary differential 
equation satisfied by the two characteristics, x = ut and y = vt, which is 

dv 6uv + v 2 + v fn 

— = — 7. . {2.6b) 

au u + 6uv + u 

The integration of ( 2.35 ) is not as simple as that of ( [2.33D since for the latter the integration 



is simply a quadrature whereas for the former a first order ordinary differential equation 
has to be solved. (This is more or less the same as the situation as the integration of a 
second order ordinary differential equation with the two symmetries ( 2.32j ) p8| .) In terms 
of the original variables the invariant is 



J 



t(x - y) 2 + {x + y) 
1 

{xy)2 



(2.36) 



in which we see that the coefficient of the independent variable, t, is simply J 2 . In principle 
we could solve ( |2.34j ) and fl2.36|) for x and y as functions of t. 

The Lie point symmetries of the system (|2.2|) are the equivalent of generalised symme- 



tries at the second order level and consequently there will be an infinity of them. We have 
no intention of presenting them all here! We do note, however, that, if G is a Lie point 

This is the equivalent result to that found 
. We can limit the number 



symmetry of (2.2), then so also is f{I,J)G. 
in the case of second order ordinary differential equations 
of symmetries found by simply making an ansatz for the functional dependence of the 
coefficient functions. Thus, if we assume at most quadratic dependence on the variables, 
we obtain the three symmetries 

d_ 

dt 

d d d 

(2.37) 



G,=- 



G 2 



G, 



dt X dx dy 
d 



+ 6xy) — + (6xy + y 2 



and. 



d_ 

dx ' v a ' " ' dy 
if we allow cubic dependence, the further four symmetries 



Ga 



Gfi 



G' 7 



t 



d_ 
at 
~d_ 
Jt 

o_ 

Of 







+ (x 2 + 6xy) — + (6xy + y 



+ (x 2 + 6xy) — + (6xy + y 



+ (x 2 + 6xy) — + (6xy + y 



yf 



dx 



d_ 

dy 
d_ 
dy 
d_ 

dy 



(2.38) 



are added. We note that the power of the independent variable, t, lags behind the powers 
of the dependent variables. 
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Table 1. The four sets of values for the coefficients of the leading terms of the Right Painleve 



Series for (3.1) 



a 


P 


7 


-1 


1 


1 


1 


1 


-1 


1 


-1 


1 


-1 


-1 


-1 



3 The simplest three-dimensional system 
3.1 Painleve analysis 

From a casual observation of the three three-dimensional systems listed in the introduction 
one cannot, a priori, judge which of them is the simplest. We shall presume that it is 
( |1.2| ) on the basis that the system contains the fewest number of terms of the three. Our 
first step is to write the system in the simplest form by rescaling the variables P and Q 
so that a becomes unity. Thus we examine the system 



w 



vw 



wu 



uv 



(3.1a) 
(3.1b) 
(3.1c) 



which is a special case of the Rikitake System. 



Because of the structure of the system (3.1) there can be no subselfsimilar symmetry 
and the leading order behaviour is that of a simple pole. To make up for this simplicity 
there are four possible sets of values for the coefficients of the leading terms. If we take the 
leading terms to be qt" 1 , (3t~ 1 and 7T -1 , the possible combinations are given in Table 
[l]. The resonances are determined from the solution of the characteristic equation 

1 -7 ~P 

0. (3.2) 



"7 



"7 
1 



-a 



-a 
1 



They are r = —1, 2(2). The repetition of the resonance at r = 2 is not a cause for concern 
since the matrix for the determination of the resonances is a real symmetric matrix and 
so the geometric multiplicity of the eigen vectors equals the algebraic multiplicity of the 
eigenvalues. The expansion up to the resonance is 



(3.3) 



and so the system ( |3.1[ ) possesses the (strong) Painleve property. 
3.2 First reduction and symmetry analysis 

Differentiation of ( |3.1| a) and the substitution of ( |3.1[ b) and Q3.1| c), division by u and the 
same process of differentiation and substitution gives the third order equation 



u 




a 


-H 




" 7 " 




' P~ 




V 




P 




1 







■} 


w 




. 7 . 






_ _ 




1 





u u 



iiu — 4tm 3 = 



(3.4) 
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which, at the level of contact symmetries, only has the two symmetries relating to time 
translation and rescaling. The Painleve analysis of ( |3.4[ ) reveals three patterns of leading 
order behaviour. The all terms dominant case has a simple pole and the resonances are 
r = —1,2 and 4. There is no inconsistency at the second resonance and this case passes 
the Painleve test. Subselfsimilar behaviour can be found with the first two terms which 
share the two homogeneity symmetries, td/dt and ud/du. The leading order terms are 
in r° and r . For both we obtain a Taylor series containing three arbitrary constants. 



Consequently (^4) possesses the (strong) Painleve property. 

We reduce the order of (3.4) using u 2 and u 2 , the two invariants of G± = d/dt, as 



the new independent (x) and dependent (y) variables to obtain the second order ordinary 
differential equation 

y" = 2- (3.5) 
Equation fl3,5|) is a linear second order equation of very simple appearance. This indicates 



that (3A) has nonlocal symmetries in sufficient supply to give the necessary right point 



symmetries for the second order equation. Equation ( |3.5[ ) has the solution 

y (x) = A + Bx + x 2 (3.6) 
and so the solution of (13.41) is 



±(t-t )= f ; ^ (3.7) 

V ' J + Bu 2 + u 4 V ^ 

in which the integral can be evaluated in terms of elliptic integrals, the expression of which 
is not informative and so we omit it. Because of the cyclic symmetry of ( |3.1| ) v and w are 
given by expressions of the same form as fl3.7j ). 

It is a trivial matter to obtain the two autonomous integrals of the system ([O]) . They 

are 

I = u 2 - V 2 and J = u 2 - w 2 . (3.8) 

The third cyclic expression is not independent of I and J. The determination of the 
necessary invariant to complete the set (with I and J) using the selfsimilar symmetry, 
G2, seems to be impossible. The one invariant which comes moderately easy is simply the 
ratio of / and J which is of no use. 



As in the case of the two-dimensional system the system (3.1) possesses an infinite 
number of Lie point symmetries and we confine our attention to those symmetries which 
are up to cubic in the variables for the coefficient functions. Thus we have 
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G 5 = u 
G 6 = v 
Gy = w 
G 8 



d d d d 

— + vw— + wu— + uv— 
dt du ov dw 

d d d d 

— + vw— + wu— + uv— 
at ou ov ow 

d d d d 

— + vw— + wu— + uv— 
at ou ov ow 

2 „.2\ d 



(3.9) 



Gg - 

G\o 
Gn 



r ~ v ~ ] Ft 

(2 2\ ° 



{u 2 - v 2 ) 
{u 2 -w 2 ) 



t® + ^ + ^ + ^ 

dt du dv dw 

d d d 

dt du dv dw 



We observe that the two first integrals, I and J, occur as coefficients in Gs to Gn. 
Although we noted the possibility of this occurrence for the system ( |2.2D , the form of 
the first integral and invariant for that system was such that symmetries of that type 
would not be found using the ansatz we have adopted. In the case of the system ( |3.1| ) the 
two first integrals fitted in with the ansatz and so we see the appearance of symmetries 
containing them as common multipliers. Given the fact that we cannot determine the 
invariant for this system, the likelihood of us making an ansatz for the structure of the 
symmetry which would include symmetries containing the invariant as a multiplier is most 
unlikely. Indeed, a priori the likelihood of an Ansatz producing such a multiplier would 
be so low that its occurrence would verge on the miraculous! A point which may be worth 
future consideration is whether the symmetries associated with a particular first integral 
have the integrals/invariants as multipliers or is there simply no ordained structure as we 
proposed for the system (|2.2j)? please consider!) We recall that for the system (|2.2|) we 



did not have / but 1 2 as the coefficient of t, whereas here we have I and J themselves. 



4 The second three-dimensional system 
4.1 Painleve analysis 



In (L4) the variables can be rescaled to give the system 



u = vw 

v = v(au + v + w) 
w = —wibu + v + w), 

where a = 1 — fi/u and b = 1 — u / u. We assume that a ^ b. 
In (4.1) we set the variables at 



(4.1a) 
(4.1b) 
(4.1c) 



u = ar 



v = (3r q , w = 7T r 



(4.2) 
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and obtain the following sets of powers of r 
p- 1 q+r 

q- 1 q+p 2q q+r . 
r- 1 r+p q+r 2r 

Evidently one set of dominant behaviours given by p = q 
determined from the set of equations 

—a = /?7 

-P = (3(aa + + 7) 
_ 7 = —j(ba + P + 7) 
and we have 

2 



(4.3) 



-1. The coefficients are 



(4.4) 



a 



P + l 



a + b 



Pi 



a — a — b 1 11 a — b 

with particularly unpleasant expressions for P and 7 which we do not bother to write 
There are two other sets of possible dominant behaviours given by 



(4.5) 



-1, r>-l, 



P 



r or 



■1, 



> -1, 



V 



(4.6) 



There is a symmetry between these two possibilities and so we just consider one, the 
second. Since the powers are required to be integral, q and p cannot be negative integers 
as the analysis requires and so we must make a full substitution. We take 



w 



E 

t=0 



a,;T 



(4.7) 



i-l 



On equating coefficients of like powers of r to zero we obtain the following set of informa- 
tion 



1 



a 01 

ai 

1 -ifean '^b 2 al 



(4.8) 



(a — 6)aoai 
\{2a — b)aoa\ 

with the rest of the coefficients following naturally. We observe that ao and a\ are arbitrary 
and these, in combination with to, provide the three arbitrary constants of integration re- 
quired to give the general solution of (4.1). Thus the subselfsimilarity dominant behaviour 
leads to a general solution and not a particular (or singular, an unfortunate word in this 
context) solution. 

It remains to look at the resonances of the all terms dominant singular behaviour. We 

set 



u 

V 

w 



ar 1 + fiT r 1 



r-l 



Pt~ x + 



VT 



(4.9) 



7T 1 + OT T X , 
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observe the following simplifications 



aa + /3 + 7 = -l 6a + /? + 7 = 1 



(4.10) 



and obtain the characteristic equation for the resonances 



r — 1 —7 — (3 
-a(3 r-p -p 
67 7 r + 7 



(4.11) 



We obtain 



-1,2,0-7. 



(4.12) 



The first two values for the resonances are fine, but the third is almost certainly going to 
be a disaster except for some very specific values of the parameters a and b. Since the 
third contains the difference of the two roots of the quadratic, even for "nice" values of a 
and b the value of r could be terrible. 



We conclude that generically (4J) does not possess the Painleve property and is not 
integrable in the sense of Painleve. There is an amusing contrast with the results of the 
singularity analysis of the equations for the Mixmaster universe. In that analysis the all 
terms dominant behaviour was better than the not all terms dominant behaviour 0, 10]. 



4.2 First reduction and symmetry analysis 



If we define r = vw, in ( |4.l| a) and ( 4.17 ), we reduce (4JJ to the following two-dimensional 
system 



u = r 

r = (a — b)ru. 



(4.13) 



When we make the Painleve analysis of ( 4.131 ), we h n d that the leading order behaviour 



is 



(a-b) 



(a-b) 



(4.14) 



The resonances occur at r = —1,2 and, as one would expect from a system of such 
symmetry, there is no incompatibility at the second resonance. We obtain 



(a-b) 



T 1 + «T 



(4.15) 



(a-b) 



t * + a. 
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The Lie point symmetries of ( 4.13 ) up to the third order in the variables are 
d 



Gi 
G 2 
G 3 
G 4 

G 5 
G% 
G 7 = 21 



dt 

^ <9 d d 

dt du dr 

9 , ^ 9 
r— h [a — b)ur 



t 



du 
d 



Or 







dt +r d-u + {a - h)ur dr 
di +r du- + {a - b)Ur dr 



(4.16) 



9 9 < t x 

dt + r d~u + ( ° - 

d_ 

dt 



(21 - (a - b)u 2 ) t^+u ((a - b)u 2 



A \ 9 A 2 9 

4r) — 4r — . 

ou or 

with the simpler system ( |4.13 ) it is possible to make the computations at least to order 
four in the variables, but there is little point in writing them down as most of them are of 
the type of G4 except with more elaborate coefficients. 



4.3 Symmetries and solution 

The combination ( |4.1[ b) w + v (|4.1| c) gives 

(vw)' = (a — b)u(vw) 
which, when combined with (|4.1|a), becomes 

il = (a — b)uu. 
We note that fl4.18| ) has only the two symmetries 



c - 9 

Gl -dl 







(4.17) 
(4.18) 

(4.19) 



2 dt~^~du 



We observe that to integrate the system ( |4.1[) it is only necessary to go to a second 
order ordinary differential equation. Equation Q4.18; ) is integrated once to give 



u = I+\(a- b)u 2 

which immediately leads to the quadrature 

du 



(4.20) 



t-t 



I+Ua-b)u 2 



1 



1 



, , 1 a-b\2 ( a-b\ 2 

2(0 — b)(t — to) = I ^r-^ I arctan | —r-r- I u. 



21 



21 



(4.21) 
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inversion of ( [4.21 ) and the trivial quadrature of ( 4.17 ) immediately give 



u = AtanQ(t — to) 
vw = B sec 2 f2(t — to) 



(4.22) 



l I 

in which A = [21 /(a — b)] 2, Q = [^/(a ~ b)] 2 , B is another constant of integration and 
we have assumed that the parameters in Q4.20| ) are positive. We consider this case only 
as the other possibilities require a similar discussion and this would be repetitive. 
We are now left with the integration of (fO|c) which now has the simple form 



w + bA tan Q{t - t )w + w 2 + B sec 2 fi(t - t Q ) = 0. 



(4.23) 



Equation ( gjg) is a Riccati equation and we transform it to a linear second order equation 
by means of the transformation w = rj/rj to obtain 



rj + bA tan Q(t - t )f] + B sec 2 Q(t - t )rj = 0. 
Equation (4.24) can be improved in appearance by the change of variables 
r](t) = y{x) x(t) = tanO(i — to)- 



We now have 



(1 + 



x I y 



2b 



xy + 



2B 



I (a - b) 



y = o. 



(4.24) 



(4.25) 



(4.26) 



Equation (fO|), being a linear second order differential equation, has eight symmetries. 
This knowledge is not of much use in the solution of the equation since seven of these 
symmetries require a knowledge of the solution. The eighth is the homogeneity symmetry, 
yd/dy, which was used to obtain this equation from the Riccati equation. Equation (4.26) 
is a little disturbing because it contains two constants of integration and one would prefer 
to discuss the solution with only the parameters, a and b, present. The constants of 
integration can be removed by dividing ( [4.26| ) by y and differentiating with respect to 
x. Unfortunately the resulting third order equation looks too awful to contemplate and 
that line of action is not pursued here. In general the solution of ( 4,26j ) is in terms of 
the hypergeometric function which is a bit too diffuse for our purposes. There are certain 
values of the constants for which solutions are known ||29|] for these have been used in the 



solution of the problem of the Tikekar superdense stars [35, £3g, 37]. Kamke [20, 2.261, p 



470] provides a method for obtaining polynomial solutions, again for certain values of the 
constants. 

From the solutions of ( [4.26 ) we can find w and hence v from ( f4.22| ). Consequently the 
passing of the Painleve property depends critically on the solution of Q4.26 ). On the other 
hand there is no constraint of this type on u. 

The presence of the two parameters a and b in the system ( |4.1| ) makes the analysis of 
the Lie point symmetries of the system problematic. The symmetries up to the quadratic 
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symmetries are 
c - d 

Gl ~¥t 

_ d d d d 

^2 = -i^T + U— + V— + tt> — 

or au ov ow /a 97^ 

d d d ' 

G3 = vw— — h v(au + v + w)— w(bu + v + w)t— - 

ait aw au; 

G 4 = [(a - 6)u 2 - 2vu;] J^. 
In G4 we recognise that the coefficient of a/at is simply the first integral in ( 4.20| ). 



5 The third three-dimensional system 

5.1 Painleve analysis 

We write (|1.3[) in the equivalent form 



ii = (2u + Xv)w (5.1) 
v = (-An + 2v)w (5.2) 
w = u 2 + v 2 + w 2 . (5.3) 



We consider the Painleve analysis of the system (|5.l| ) - (|5.3| ) . We find that the powers 
in the different terms are 

— 1 r —p + q + r 

— 1 — g+p + r r (5-4) 

— 1 2p — r 2q — r r, 

where we have assumed the leading terms to be 



v = (3r q (5.5) 

W = ^T r . 

The only possible singular behaviour is p = q = r = —1. The coefficients a, j3 and 7 are 
determined from the system 

-a = (2a + A/3) 7 (5.6) 

-0 = (_Aq + 2/?) 7 (5.7) 

- 7 = a 2 + /3 2 +7 2 . (5.8) 



The combination of a(|5.6|) + /?( [5. 7| ) gives 

- (a 2 + /3 2 ) = 2 (a 2 + /? 2 ) 7 (5.9) 
from which it follows that either 7 = — i or a 2 + /3 2 = 0. We consider each in turn. 
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From 



a 2 + /3 ; 



it follows that 
2 _ 1 

4- 



With this value of 7 (5.6) and ( |5.7| ) become 
A/3 = Aa = => A = 



(5.10) 



(5.11) 



since the assumption of leading order behaviour requires that a and f3 be nonzero. A 
similar analysis for the second possibility leads to the result that A = ±i. 
Case: A = The system is now 



it = 2uw 
v = 2vw 

2,2, 2 
W = U + V + w 

and to determine the resonances we make the substitution 

U = OtT~ l + /XT r_1 

v = /5r _1 + h>T r ~ l 
w 

The characteristic equation for r is 



(5.12) 



(5.13) 



7T 1 + (JT r 1 . 



" - 1 - 27 

r- 
-2a 

which gives the values 
r = -l,0,l 





1 - 27 
-2/3 r 



-2a 
-2(3 
1 - 27 



(5.14) 



(5.15) 



when the values of the constants are substituted. There is only one eigenvector which is 

(5.16) 




That there is only one eigenvector indicates that it is necessary to introduce a logarithmic 
term and the for the entry of w indicates that the logarithmic term must come in with 
the other two variables, u and v. 

Case: A = ±i With the same substitution as that in ( |5.13j ) the characteristic equation 
for the resonances is given by 



r - 1 - 27 
A7 r 
-2a 



-A7 2a + A/3 
I-27 -Aa + 2/3 
-2/3 r - 1 - 27 







and with the constraints on the constants gives the values 
r = -2,-1,0. 



(5.17) 



(5.18) 
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The value r = is the critical one. Two possible results occur. In the case that a = \(3 
fj, and v are arbitrary because the rank of the matrix becomes one. There is no need to 
introduce a logarithmic term. In the case that a 7^ A/3 this is not the case because \x and 
v are related by 2afi + 2{3v = 0. A logarithmic term must be introduced. 

We conclude that only for the specific values of A shown above does the system (5.1) - 
(|5.3|) pass the Painleve test. 



5.2 The two reductions 

We now consider the integration of the system ( |5,lD - ( |5.3| ). We observe that the combi- 



nation u (|5.l| ) + f (|5.2|) gives 



uii + vv = 2(u 2 + v 2 )w (5.19) 



so that differentiation of ( |5.3| ) with respect to time gives the second order ordinary differ- 
ential equation 

id - 6ww + 4u7 3 = (5.20) 

which we recognise as the equation we already met in ( |2.10| ). We know that it possesses 
eight Lie point symmetries, that under a Riccati transformation it becomes the simplest 
third order equation and that the solution is 

B + Ct . 
W = -A + 2Bt + Ct*- (5 ' 21) 



Now that we know the functional form of w(t) we see that ( |5.1[ ) and (|5.2| ) now constitute 
a linear system. This system is rendered autonomous by the introduction of a new time 
variable ||, ||] defined by 

= w & T = log (A + 2Bt + Ct 2 ) ~ 2 . (5.22) 

at 

If we denote differentiation with respect to r by a ', the linear system is 



u \ _ ( 2 A \ / u 
v ~[ -X 2 [ v 



(5.23) 



The solution of ( 5.22j ) is 



U \ I -l \ i( 2 +i A)r , / M „p( 2 - iA ) 



; y J . ^ z+lA > T + 1 J ve (Z -^> T . (5.24) 

There appear to be too many constants of integration, but the three in the expression 
for w are in reality two and a small calculation shows that the constants of integration in 
( 5.24| ), viz /x and u, are related to the other constants via 



2[iv = B 2 -AC (5.25) 



so that the number of arbitrary constants is in fact three. In terms of the variable r the 
solutions for u and v are analytic. Given the definition of r in (5.22), it is evident that 
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the solution in terms of the original variable, t, is not analytic in general. In terms of the 
variable t the solution for w has simple poles as the only singularities. If we express this 
solution in terms of the new time, we find that 



w 



±e 2r \/ .B 2 — AC + Ce — 2t 



(5.26) 



which has a branch point singularity. Of course, to define a new time in terms of the 
solution of one of the variables means that the problem has been reduced from a three- 
dimensional system to a two-dimensional system and a linear one at that. 

In fact we can interpose an intermediate system between ([h^) and ( |5.20| ) by introducing 
a new variable, r, defined by 



2 2,2 

r = u + v . 



The system is 
r = 2rw 

2 2 

w = r + w . 
The solution to the system ( 5.28Q is ( |5.21 ) and 
B 2 -AC 



(5.27) 



(5.28) 



(A + 2Bt + ct 2 y 



(5.29) 



in which the additional constant of integration has already been identified in terms of the 
three constants found in the solution for the function w. 

We now consider the Painleve analysis of ( [5.20D , ( 5.28 ) and fll.3j ) (in the present vari- 
ables) in turn. In the case of ( |5.20| ) the only possible singular behaviour is given by p = —1, 
but there are two possible families of solutions since the coefficient of the leading term can 
be either — 1 or — i. To determine the resonances we write 



w = ar 1 + /3r r 1 
so that we can treat both cases simultaneously. We find that 

r = -1,6a + 4 
so that we have 



a 



i 



■1,1 



o 



-1 



-1,-2. 



(5.30) 



(5.31) 



(5.32) 



Evidently the first gives the right Painleve series, that is the Laurent expansion in the 
neighbourhood of the singularity, and the second the left Painleve series, the asymptotic 
expansion, that is the Laurent expansion away from the singularity. 

For the system ( |5.28| ) the solution for w is given in (5.21) and it is a simple quadrature 
to find that 



B 2 - AC 
A + 2Bt + Ct 2 ' 



(5.33) 
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We determine the leading order behaviour by writing 

r = ar p w = (3r q (5.34) 
and substituting this into (5.28). We have 



a P T p 1 = 2a[3rp + q 
Pgr*' 1 = a 2 T2P + f3 2 T 2q . 

Evidently there are two possible behaviours. We can have p = q = —1 or p > —l,q = —l. 
The roles of p and q cannot be reversed in the second possibility. We find that a = ±| and 
/3 = — \. The resonances are given by r = ±1. The first two terms of the right Painleve 
series are 

where fi is the second arbitrary constant required for the general solution. For the second 
possibility we set 

r = ^ ay w = Y, biT 1 - 1 (5.37) 

i=0 i=0 

since the standard Painleve analysis breaks down in the case that the leading order be- 
haviour is not singular in all variables 0, . We substitute (|5.371 ) into ( |5.28| ) to obtain 

E^-^EE^^" 1 

(5.38) 

5> - m^ 2 = E E + E E hb^-\ 

i=Q i=0 j=0 i=0 j=0 

When we compare coefficients of like powers of r, we find that all coefficients except bo 
vanish so that the solution is a singular (particular) solution and is unrelated to the general 
solution. The same result is found if we commence the series for r at a higher power of r. 
Thus we see that the Painleve property holds equally well for ( 5.28|) as it did for ( |5.20 ). 



Bearing in mind the results of the Painleve analysis of the system (|ll|) we see that the 
decomposition of the system (5.28) into the system (1.3) and even the possibility of the 



preservation of the Painleve property is restricted to those special cases for the value of A 
obtained above. 

In the case of this system going from a system of two equations to a system of three 
equations has had a drastic effect upon the possibility of the possession of the Painleve 
property. Generically, it for arbitrary values of A, this system of three equations does not 
have the Painleve property. 



5.3 Symmetries 

The Lie point symmetries for ( [5.20 ) are 
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G t +w 

dt dw 

( 9 2 9 N 
G3 = W {dt +W d^ / 

d d 
Gi = tw— + w 2 {l + tw) — 
dt dw 



G 6 
G 7 
G s 



dt dw 

d d 
+2 -n — + w(l + 2tw + 2t 2 w 2 ) — 
dt x ' dw 

Q , o + „..\ A , fft + 2„„2 . „ + 3„„3\ 9 



„ n d 

G ^ =t W dt 

= t 2 (3 + 2tw) ^ + (6t 2 w 2 + 4t 3 w 3 ) ^- 

d 

3tw + 2t 2 w 2 ) 7— 
' dw 

i... , hj-2„..2 , oj.3„..3\ 9 



9 I 
dl + tW ^- 

d 

= t 3 (1 + tw) — + t (1 + 3tw + At 
those of ( 5.2$ ) up to cubic in the variables are 



= t 2 (i + H|- 

! 2 w 2 + 2t 3 w 3 ) 

dw 



r - 9 



t 9 + 9 + 9 

dt dr dw 



° 2 ~ ~ L di^ 

^ dw 



9 . . > , 2 
" + w 



■ + (r £ 



d 



G 4 = t 
G 5 = r 
G e = w 



d 



G 3 = Arw— 
Or 

— + Arw — + (r 2 + w 2 ) — 
dt dr dw 

9 A 9 12 2\ 9 

+ Arw— + (r + w ) — 
or Ow 



dt 

'8,8, 
— + irw— + (r 
dt dr 



d 



and of the original system, ( |1.3| ), (in the present variables) 

r - 9 



are 



d d 



d d 

G 2 = -t— + u— + v— + w— 

dt du dv dw 

d d d 

G 3 = (2u + Xv)w— + (-Ait + 2v)w— + (u 2 + v 2 + w 2 ) — 
du dv dw 

d d d d 

= t — + (2u + Xv)w— + (-An + 2v)w— + (u 2 + v 2 + w 2 ) — 
dt du dv dw 

' d d d 

— + {2u + Xv)w— + (-Xu + 2v)w— + (n 5 
at du dv 

d d (9 

— + (2u + Xv)w— + {-Xu + 2v)w— + {w 
dt du dv 

^- + (2u + Xv)w-^- + (-Xu + 2v)w^- 
dt du dv 



G 4 = t 
G 5 = u 
G 6 = v 
Gj = w 



L 2 + V 2 +W 2 )^- 

dw 

2 2\ 9 ' 

+ v 2 + w 2 ) — 
dw 

d n 



9 i 2 2 2\ 9 

-+(n 2 + , 2 + . 2 )- 



(5.39) 



(5.40) 



(5.41) 
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The symmetries in ( [5,40 ) and in ( |5.4l| ) are essentially the same with the latter having an 



additional symmetry due to the presence of the additional variable. 

6 The four-dimensional system 
6.1 Painleve analysis 

The parameter /3 in ( |1.1| ) is an essential parameter which cannot be removed by rescaling. 
In fact any attempt to remove it will simply lead to the appearance of another parameter 
in the other terms. We rewrite the system as 

it = u(u — w) — vx (6.1a) 

v = v[bu + (b - 2)w] (6.1b) 

w = —w(u — w) — v x (6.1c) 

x = -x[(b-l)u+ (b-3)w]. (6.1d) 

We perform the Painleve analysis of (|6.1|). In the case of all terms being dominant the 
singularity is a simple pole and the coefficients of the leading order terms must satisfy the 
following set of equations 

—a = a 2 — a-y — 

-/3 = 6a/3 + (6-2)/3 7 

—7 = — a 7 + 7 — j3o 

S = -(b-l)aS-(b-3)jS. 

Since the leading order behaviour assumes that the coefficients are nonzero, the second 
and fourth of fl6.2|) give 



6a+(6-2) 7 = -l 
(6 -l)a + (6-3)7 = 1 

from which it follows that either a = 7 = — 1 or 

q + 7 = — 2 and + 7 = — 1 (6-4) 

which is a contradiction. Substitution for the value of a and 7 into ( |6.2| ) requires that 
the product /3 7 = — 1 and, more importantly, b = — 1. Thus the Painleve analysis of ( |6.l| 



cannot even begin without a restriction of the value of b. We observe in ( 6.25| ) and ( |6.26 ) 



that this value of b removes the possible singularity in the trigonometric term, but we are 
still left with the square root to give the solution a branch point singularity. 

Nevertheless for the sake of completeness we analyse the situation for this particular 
value of b. On making the usual substitution for the leading order behaviour we obtain 
the following set of exponents 



(6.5) 



p- 


1 


2p 


p + r 


q + s 


q - 
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p + q 


q + r 




r — 
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p + r 


2r 


q + s 


s — 


1 


p + s 


r + s 




-1 




P 
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in which the separated bottom line gives the essence of the content of the second and 
fourth lines. When all terms are dominant, the singularity is of order —1. However, there 
are various options for subdominant behaviour. These are listed below as 



p q r s 

-1 -1 -1 -1 

-1 > > > 

-1-1 >0 >0 

-1 > > -1 



Comment 
generic 

plus 
variants 
p <-> r. 



(6.6) 



We note that the leading order behaviour in the subdominant cases must contain either p 
or r. 

Only in the case of all terms dominant can we apply the standard algorithm. The 
coefficients of the leading order terms are found from the solution of the system 



—a = a — «7 — (35 
-2(3 = 3a(3 - (35 
—7 = cry + 7 2 — (35 
-25 = -a5 + 3-f5 



(6.7) 



and are 



a 



7 



-1 (35 



ie one of the coefficients is arbitrary, 
solution of the system of equations 



( r 5 

-3(3 2r 

-1 7 

\ 5 



-1 \ 

(3 

r (3 

-35 2r J 



a 

\ P J 



(6.8) 

The resonance behaviour is determined by the 






V / 



(6.9) 



The resonances are given by r = —1(2), 0,2. The double resonance at —1 means that 
we cannot have the general solution since there will be at most only three arbitrary con- 
stants. For the resonances at r = and r = 2 we find that the corresponding vectors are 
respectively 



ao 



"I \ 
±1 

1 

v I 1 / 



1 1 \ 



and 3l2 



(6.10) 



V + 3 7 



Thus we see that the resonance r = indicates the existence of two possible solutions. Each 
of these solutions will contain an arbitrary constant introduced at the r = 2 resonance. 

In the cases of not all terms being dominant the standard algorithm cannot be applied 
and so we must substitute Anzatse appropriate to the particular pattern desired. For the 
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second set of indices in (|6.6| ) we write 

E^ 1 

8=0 

E 6 * 



i=0 



r - > OjT 

^ . (6.H) 

i=o 



8=0 



«0 


= -1 


h 


= 


CO 


= 


f/ 


= 


a I 


= 


h 


= 


Cl 


= Cl 


dl 


= 


«2 


= lc 


h 


= 


C2 


= 




= 


"3 


= 


b 3 


= 


C3 


= 3" C 1 


d 3 


= 0. 



and after equating the coefficients of like powers of r to zero we obtain the following set 
of coefficients 



(6.12) 



Evidently the functions v and x are identically zero for this pattern of singularity. 

For the third and fourth sets of indices in (|6.6| ) we obtain the same expansion as in 
( |6.12j ). All the solutions for the different possible singularity structures are particular 
solutions and, as we expected above, the system fl6.1|) does not satisfy the requirements of 
the Painleve test for any of the values of the parameter b. 

6.2 The general solutions 

From the combinations of ( |6.1| a) — ( |6.1| c), fl6.1f j)x + v ( |6.1| d) and (|6.1|a)w + u (|6.1b ) we 
obtain respectively 

u — w = u 2 — w 2 (6.13) 

(vx)' = (u + w)vx (6-14) 

(uw)' = —{u + w)vx. (6.15) 



From the combination of ( 6,14| ) and ( 5.15| ) we obtain the first integral 

—I 2 = uw + vx (6.16) 



and of ( |6.13D and ( |6.15 ) the first integral 



J=(u — w)vx, (6-17) 



where I and J are constants of integration. In ( |6.16 ) we have chosen the label for the first 
integral as —I 2 for later convenience. The other cases which we do not treat here can be 
treated in a fashion similar to what we are about to do. 
When we substitute ( 6.16| ) into ( |6.1| a), we obtain 



u 2 + I 2 (6.18) 



Symmetry, Singularities and Integrability in Complex Dynamics I 



471 



which can be integrated by an elementary quadrature to give 



u(t) = Jtan(i — to), 



(6.19) 



where to is the constant of integration. By the elimination of vx between ( |6.16| ) and ( |6. 1T|) 
we obtain a quadratic equation for w(t) which, when we substitute for u(t) using ( |6.19 ), 



gives 



w(t) 



sin2(t — to) 



{-Jcos2(t - t ) ± \J[P + J sin 2 (t - t )(l + cos(t - t ))/I}} 

(6.20) 



so that ( 6. 17| ) gives 

-±/sec 2 (t - t ) {/ ± ^[P + Jsin2(t - t )(I + cos(t - t ))/i]} 



vx 



(6.21) 



We observe that the parameter, b, of the system ( |6.1[) does not appear in the solutions 
obtained above. The presence of the square root suggests that contrary to the comment 
of Golubchik and Sokolov that the system (6.1) probably does not pass the Painleve- 
Kowalevskaya test for generic values of the parameter is something of an understatement. 
To determine the functions v(t) and x(t) we can now treat (p\T|b) and (|6.1| d) as linear 
equations for the two variables. We obtain 



and 



v(t) = ETexp 



x(t) = K exp 



(bu + (b - 2)w))dt 



((6- l)it + (b - 3)w))dt 



from which it is evident that 



vx(t) = KK exp 



'u + w)dt 



(6.22) 



(6.23) 



(6.24) 



where the notation vx(t) is meant to indicate vx as a function of time, which gives us a 
neat way to write down the solution without going to the effort of evaluating the integral 
in either ( 6.22j ) or ( |6.23 ). Given the appearance of the square root in w one would not 
expect to be able to evaluate either integral with ease. Using this device we find that 



and 



v(t) = K 



x(t) = K 



vx{t) 
KK 



vx{t) 
KK 



b~2 



sec 2I (t - t ) 



3-6 



cos 2/ (t - t Q ). 



(6.25) 



(6.26) 



The exponents in the expressions in ( 6.25| ) and (6.26) definitely indicate that one should 
not expect the system, fl6,ip, to satisfy the Painleve test for generic values of the parameter, 
b. 
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6.3 Second reduction 

Another way to reduce the order of the system Q6.1| ) is by the following procedure. We 
define the variables 

s = u — w 

y = u + w (6.27) 

z = vx 



which, with the use of the system ( |6.1[) , satisfy the three-dimensional system of first order 
equations 

s = sy 

y = s 2 -2z (6.28) 
z = yz. 

Consider the system fl6.28|) . When we make the usual substitution for the leading order 
behaviour, we require the following terms to balance: 

aprP- 1 = af3r p+q 

pqri- 1 = a 2 T 2p - 2 7 r r (6.29) 
jrrr — 1 = (3jT q+r . 

From the first and last of ( |6.29| ) it is obvious that q = — 1. If all terms are to be dominant 
we also have p = — 1 and r = —2. However, when we solve for the coefficients, we find that 
the first of ( |6.29| ) gives (3 = — 1 and the third (3 = —2 unless both a and 7 are zero. This 
immediately implies that f3 = which contradicts the idea of leading order behaviour. 
Consequently we must conclude that the case of all terms dominant is impossible. 

There are two possible sets of subdominant behaviour. These are p = —l,r > —2 and 
p > —1, r = —2. Effectively these two cases are p = —1, r = —1 and p = 0, r = —2. 

We commence with the second and, since one of the putative solutions is not singular, 
we make the ansatz 

■s - > a 7 -T J_1 



(6.30) 



i=0 
i=0 



E 



CiT 





and substitute these into ( 6.28 ). After a modicum of calculation we find that 



«0 


= ±1 


b 


= -1 


CO = 


h 


a 1 


= ±b x 


h 


= h 


Cl = 


m 


a 2 


- ± 2 -b 2 




_ 5°1 


C2 = 


-I 


«3 


= ±|6? 


h 
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(6.31) 
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in which it is quite obvious that there are only two arbitrary constants, the location of 
the singularity and the coefficient b\. Consequently the Laurent series does not represent 
the full solution in the neighbourhood of the singularity, but a particular solution of the 



system ( 6,28 ) 



For the first possibility of subdominant behaviour we substitute 



s = 


\ — ^ 

E 

i=0 


diT l 








y = 


i=0 


x. i— 1 








z = 


E 

i=0 










and find 
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CO 


= -1 


ai 
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Cl 
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-2c 2 


C2 


= c 2 
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C3 
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a. i 
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5 C 2 


C4 


_ 3 
5 



(6.32) 



(6.33) 



with only the even coefficients and C2i being nonzero. The solution depends upon 
two arbitrary constants, but is definitely very particular for the system since one of the 
functions is identically zero. It is probably incorrect to regard this solution as the general 
solution for the two functions y(t) and z(t). 

Somewhat dismayed by the result displayed in ( 6.32| ) we are smitten by a stroke of 



genius and think to apply a left Painleve series to see if we can obtain a solution which 
does contain the function s(t). Thus we write 



-i-2 
OiT 



E' 

i=0 

2 biT '^ 1 (6.34) 



and find 



i=0 

i-2 



CiT 

1=0 



«0 


= ±(26 2 - 


6f)3 


6o 


= -2 


c = -1 
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bi 


= 61 
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= ±i(26 2 


-&?)*(&? -fig) 


6 2 
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C2 = 3 (^2 - 


«3 


= ±g(26 2 


-6?)5(f6i62-26?) 


63 


= ^^1^2 - Jft? 


C3 = -56162 



(6.35) 



in which both b\ and 6 2 are arbitrary. However, this is a Laurent expansion at infinity 
and so the expansion is in negative powers of the variable t and so we have a particular 
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asymptotic solution and not the general solution of ( 6.28j ). We recall that the single third- 



order equation ( 6.37| ) had partial expansions for both left and right series for one of the 



values of the coefficient of the leading order power. Nevertheless we are inclined to think 
that this is possibly the first time that a left Painleve series indicating a particular solution 
has been demonstrated in the case of subdominant behaviour. 

By differentiating the second of fl6.28| ) and utilising the other two equations we obtain 



y-yy = S y. (6.36) 

With a repetition of the same procedure on ( |6.36| ) we obtain a third order equation for 
the variable y, viz 

y y -yy - 3y 2 y + 2y 3 y = o. (6.37) 



Equation ( 6.37 ) possesses just the two symmetries due to invariance under time translation 
and rescaling, even if the computation is carried out for contact symmetries. If we perform 
the Painleve analysis on ( |6.37| ), we find that the only admissible singular behaviour is a 
simple pole. There are two values for the coefficient of this term, viz a = —1,-2. For 
the first of these the resonances are given by r = — 1, 1, 2 and so, in conjunction with the 
two symmetries, ( |6.37| ) passes the Painleve test. In the case of the second value of a the 
resonances are r = —1, ±2. Consequently there is no possibility of the Painleve test being 
satisfied for a = —2. If we take the value +2, we obtain a right Painleve series containing 
only two arbitrary constants. If we take the value —2, we obtain a left Painleve series also 
containing only two arbitrary constants. One of these would be a particular solution in 
the vicinity of the singularity and the other a particular asymptotic solution, both in the 
sense defined by Cotsakis and Leach |Tl[|). 

This immediately raises the question of the integrability of ( |6.37[ ) since it does not 
pass the Painleve test for all possible singularity patterns. If we divide ( |6.37D by y 2 , the 
equation becomes exact and we obtain the first integral 

-L 2 = ^-3y + y 2 , (6.38) 

y 

where we have taken the expression for the left-hand side to be more convenient for later 
usage. We can rewrite ( |6.38| ) as the differential equation 

y-3yy + y 3 + L 2 y = 0. (6.39) 

We observe that the constant of integration, L, can be scaled out of ( 6.39| ) without affecting 



the other terms. We also observe that ( 6.39 ) is simply a variation on the equation we met 



in §§2 and 3 (( [2.10 ) and ( |5.20| )). The presence of the linear term does not affect the 
number of symmetries, which remains at eight and so implies that there exists a point 
transformation which linearises the equation, although it does make the transformation 
look more complex. 

To solve (l6T39l ) we introduce the Riccati transformation 



y = -- (6.40) 
V 
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to obtain the third order linear equation of maximal symmetry, 

V/ + L 2 r/ = (6.41) 
with solution 

rj(t) = A + BsmL(t -to) (6.42) 
so that 

V(t) = B ^('-'°> (6 .43) 
A + B sm L(t — to) 

which, provided < 1, obviously has only poles as singularities. 



Consequently we have a counterexample to the claim by Tabor [34, p 330] that the 
possession of the Painleve property requires that all possible patterns of singular behaviour 
pass the Painleve test. 



We can make a few other observations. We have noted that ( 6.3S| ) possesses eight 



Lie point symmetries with or without the presence of the linear term. As in the case of 
(|2TTo| ) and (gjgjj) , (|6.39|) without the linear term possesses both a left Painleve series and 



a right Painleve series. This is not the case when the linear term is present. Then the 
only possibility is a right Painleve series. This does not affect the integrability of ( |6.39[) 
in terms of functions with only poles as moveable singularities. The existence of the same 
algebra of Lie point symmetries means that there exists a point transformation of ( |6.39 ) 



which removes the linear term. That there is a change in the nature of the singularity 
behaviour means that this transformation cannot be homeographic. Naturally it must 
have only polelike singularities since both solutions possess the Painleve property. 
This concludes the discussion of the fourth order system. 



7 Discussion 

From the analyses of the systems of equations ( |1 . 1| ) , ( |1.2| ), (|1.3|), ( |1.4|) and ( fLlf ) several 
conclusions can be drawn. 

The development of the systems of equations, which in the original paper were 
obtained as realisations of some kinds of nonassociative algebras, considered in this paper 
may be viewed as a process of division of a single higher order equation. In the case of 
the two-dimensional quadratic system, (1.5), one can easily imagine it as being derived 



from a second order equation by a standard process of reduction to a system of two first 
order equations. It is not surprising that there was no loss of the Painleve property in this 
splitting. In the case of (|1.2j ) the three-dimensional first order quadratic system was derived 
naturally from a third order nonlinear equation. Although this equation fl3.4| ) had only 
two contact symmetries, reduction of order brought it to a linear second order equation 
with eight symmetries and consequently ([O]) must have seven nonlocal symmetries with 
nice reduction properties. It is not surprising that fl3.l|) possesses the Painleve property. 
In the case of the equations, ( |1.3|) and (|1.4j ), there is a distinct change in behaviour since 
the parameters which are present in the equations are essential parameters and cannot be 
removed by rescaling. The basic equation is really a nonlinear second order equation and 
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one could imagine this being split into two first order equations followed by a further split 
of an arbitrary nature of one of the equations into two first order equations containing the 
arbitrary parameter. This situation is repeated in the case of at a higher order in that 
one starts from a third order nonlinear equation and obtains three first order equations 
which have the Painleve property and which do not contain a parameter. The fourth first 
order equation is obtained by a parameter-dependent split of one of these three first order 
equations and the Painleve property is lost. 

One can envisage a situation in which a prediction of the preservation of the Painleve 
property can be made. Given an nth order scalar ordinary differential equation possessing 
the Painleve property, one could expect to be able to make a system of n first order 
ordinary differential equations by some sort of natural reduction and the system would 
still possess the Painleve property. Any attempt to have more than n first order equations 
by some sort of further splitting may inevitably lead to a loss of the Painleve property. It 



is interesting that the second order equations possessing the Painleve property, viz ( [2.10 ), 
(PT2CD and Q, which can be viewed as the origins of the systems ( |1.2|) , ([O^) and (|1.5|) , 
respectively are essentially the same equation and that the third order equation from 
which the system (|1 . 1|) was derived is closely related to that equation. This does pose 
some interesting questions about the basis of the work of Golubchik and Sokolov [17] from 
the point of view of the algebraic structure of scalar ordinary differential equations of the 
second and high orders. 

A third question which has only been suggested here is the question of the algebraic 
structure of the symmetries of the first order systems. There seems to be some structure 
in the symmetries listed. We note that this structure appears to be independent of the 
possession of the Painleve property which is not surprising. Algebraic properties do not 
require an analytic basis. Thus a system can be integrable in the algebraic sense and not 
in the sense of Painleve. However, because symmetry is at the basis of the process of 
reduction, a system which is integrable in the sense of Painleve must also be integrable in 
the algebraic sense. 

Finally we mention the interesting observation that a left Painleve series can be found in 
subdominant behaviour. Related to that is the whole question of the Painleve integrability 
of fl6.37l ). The solution (|6.43| ) is manifestly possessed of only poles as singularities. The 
equation has the Painleve property for one possible variety of leading order behaviour, but 
not for the other. This makes one wonder just how much is required for one to conclude 
that there exists a solution expressible in terms of a Laurent expansion about a polelike 
singularity. 
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